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Abstract 

This paper is concerned with the Riemann problem for the isentropic Chaplygin gas magnetogas¬ 
dynamics equations and the formation of delta shocks and vacuum states as pressure and magnetic 
field vanish. Firstly, the Riemann problem of the isentropic magnetogasdynamics equations for Chap¬ 
lygin gas is solved analytically. Secondly, it is rigorously proved that, as both the pressure and the 
magnetic field vanish, the Riemann solution containing two shock waves tends to a delta shock solution 
to the transport equations, and the intermediate density between the two shocks tends to a weighted 
(5-measure which forms the delta shock; while the Riemann solution containing two rarefaction waves 
tends to a two-contact-discontinuity solution to the transport equations, the termediate state between 
the two contact discontinuities is a vacuum state. 
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1. Introduction 


In this paper, we are concerned with the system of conservation law governing the one-dimensional 
unsteady simple flow of an isentropic, inviscid and perfectly conducting compressible fluid subjected to 
a transverse magnetic field (see [10, llj): 


pt + {pu)x = 0, 

{pu)t + (p + pu^ + /‘2p)x = 0, 


( 1 . 1 ) 


where p > 0, u, p, B and p > 0 represent the density, velocity, pressure, transverse magnetic field and 
magnetic permeability, respectively; p and B are known functions defined as 
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and B = k 2 p, where ki and fc 2 are positive constants. The independent variables t and x denote time 
and space, respectively. The adiabatic exponent in (1.2) can be viewed as 7 = —1 by comparing with 
the state equation p = with 7 > 1 for the polytropic gas. The gas (1.2) whose adiabatic constant 
7 = — 1 is usually called as the Chaplygin gas. 

For the isentropic Chaplygin gas Euler equations, Brenier [1] firstly studied the 1-D Riemann 
problem and obtained solutions with concentration when initial data belong to a certain domain in 
the phase plane. Furthermore, Guo, Sheng and Zhang [ 6 ] abandoned this constrain and constructively 
obtained the global solutions to the 1-D Riemann problem, in which the (5-shock developed. Moreover, 
they also systematically studied the 2-D Riemann problem for isentropic Chaplygin gas equations. For 
the 2-D case, we can also refer to [9] in which D. Serre studied the interaction of the pressure waves 
for the 2-D isentropic irrotational Chaplygin gas and constructively proved the existence of transonic 
solutions for two cases, saddle and vortex of 2-D Riemann problem. Recently, Sheng, Wang and Yin 
[13] and Wang [15] studied the Riemann problem for the generalized Chaplygin gas and obtained the 
solutions to the Riemann problem and the interactions of elementary waves. The Riemann solutions 
to the transport equations in zero-pressure flow in gas dynamics were presented by Sheng and Zhang 
in [14], in which delta shocks and vacuum states appeared. 

In related researchs of the (5-shock waves, one very important and interesting topic is to study the 
phenomena of concentration and cavitation and the formation of (5-shock waves and vacuum states in 
solutions. In earlier paper [4], Chen and Liu [4] studied the formation of 5-shocks and vacuum states of 
the Riemann solutions to the isentropic Euler equations for poly tropic gas as e —>■ 0, in which they took 
the equation of state as P = ep for p = p'^ (7 > !)• Further, they also obtained the same results for the 
Euler equations for nonisentropic fluids in [5]. The same problem for the the isentropic Euler equations 
for isothermal case was studied by Li [7], in which he proved that when temperature drops to zero, the 
solution containing two shock waves converges to the delta shock solution to the transport equations 
and the solution containing two rarefaction waves converges to the solution involving vacuum to the 
transport equations. Then, the results were extended to the relativistic Euler equations for polytropic 
gas by Yin and Sheng [17] and for Chaplygin gas by Yin and Song [18], the isentropic Euler equations for 
the generalized Chaplygin gas by Sheng, Wang and Yin [13] and for modified Chaplygin gas by Yang and 
Wang [16], the perturbed Aw-Rascle model by Shen and Sun [12], the isentropic magnetogasdynamics 
equations for polytropic gas by Shen [11], the generalized pressureless gas dynamics model with a scaled 
pressure term by Mitrovic and Nedeljkov [ 8 ], etc. 

In this paper, we study the Riemann problem of the isentropic magnetogasdynamics equations for 
Chaplygin gas and the formation of delta shocks and vacuum states as pressure and magnetic field 
vanish. The organization of the paper is as follows. In Sections 2 and 3, the Riemann problems for 
the isentropic Chaplygin gas magnetogasdynamics equations and the transport equations are analyzed 
by characteristic analysis. In Sections 4 and 5, we investigate the formation of 5-shocks and vacuum 
states of the Riemann solutions to the isentropic magnetogasdynamics equations for Chaplygin gas as 
pressure and magnetic field vanish. 

2. Riemann problem for system (1.1)-(1.2) 
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In this section, we discuss the Riemann solutions of (1.1) and (1.2) with initial data 


{p,u){x,0) = {p±,u±), ± 0 ; > 0, 


( 2 . 1 ) 


where p± > 0 and u± are arbitrary constants. 

For smooth solution, system (1.1) is equivalent to 


u p 



'Up' / p u 



= 0 , 


( 2 . 2 ) 


where 'w = {P -\- is the magneto-acoustic speed with c = {p'{p)Y^'^ as the local sound speed 

and b = {B^{p)/ppy^^ the Alfven speed. Here, prime denotes differentiation with respect to p. The 
eigenvalues of system (1.1) and (1.2) are 


Ai = M — 


fci , fejp 
“2 ' ’ 
P P 


\2= U + 


k\ k^p 


Therefore, system (1.1) and (1.2) is strictly hyperbolic for p > 0. 
The corresponding right eigenvectors are 


it = (-P, 

By simple calculation, we get 


fei klp^T 


) , it = (p, 


fei kip T 


VXi-rt = 


Skip 




/O, 1 = 1 , 2 . 


Therefore, both the characteristic fields are genuinely nonlinear. 

Since system (1.1), (1.2) and the Riemann data (2.1) are invariant under stretching of coordinates: 
{x,t) {ax, at) {a is constant), we seek the self-similar solution 

{p,u){x,t) = {p,u){i), ^=j- 

Then Riemann problem (1.1), (1.2) and (2.1) is reduced to the following boundary value problem of 
ordinary differential equations: 


-CPe + {pu)s = 0, 

-apu), + {-^+pP + ^f)^ = 0, 

with (p,m)(±cx 3) = (p±,M±). 

For any smooth solution, system (2.3) can be written as 

u-^ p 

—i'u+^ + P -\—^ —^p+2pu 


(2.3) 



= 0 . 


(2.4) 


It provides either general solutions (constant states) 


{p,u){^) = coi^stant (p > 0) 



or singular solutions called the rarefaction waves Ri and R 2 which denote, respectively, 1-rarefaction 
waves and 2-rarefaction waves. 


Ri : 


and 


R 2 : 



(2.5) 


ds, p<p-, 


ds, p> P-. 


Differentiating the second equation of (2.5) with respect to p yields Up = — 
subsequently. 


( 2 . 6 ) 


< 0, and 


4k 1 _|_ ^2^ 
p'‘‘ M 


2 p= 


fci , fcfp 
717 H- 77 - 


> 0 , 


which mean that the 1-rarefaction wave curve Ri is monotonic decreasing and convex in the (p, u) plane 
(p > 0). Similarly, from the second equation of (2.6), we have Up > 0 and Upp < 0, which mean that 
the 2-rarefaction wave curve R 2 is monotonic increasing and concave in the (p, u) plane (p > 0). In 
addition, it can be verified that lim u = -l-cxa for the 1-rarefaction wave curve Ri , which implies that 

p —^ 0 ”^ 

Ri has the -u-axis as its asymptotic line. It can also be proved that lim u = H-oo for the 2-rarefaction 

p—^+00 

wave curve i? 2 . 

For a bounded discontinuity at ^ = cr, the Rankine-Hugoniot conditions hold: 


-a[p] + [pu] = 0, 

-a[pu] + [-^-f+pu^+^^] = Q, 
where [p] = p — p-, etc. Solving (2.7), we obtain two shock waves Si and S 2 


k-i 



)(P-P-), P> P-, 


)(P-P-), P<P-- 

Differentiating the second equation of (2.8) with respect to p, for p > p- we have 


1 


y pp- pp- 


+ 


2 p 


/ 2A:i ^2 P2 
Vp-p3 p-U 2pu 


ko ^ ^2p— 


kRp+p.)\ P-P 2pp 2p2p 


< 0 , 


(2.7) 


( 2 . 8 ) 


(2.9) 


which mean that the 1-shock curve Si is monotonic decreasing in the (p, u) plane (p > 0). Similarly, 
from the second equation of (2.9), for p < p_ we have Up > 0, which mean that the 2-shock wave curve 
S 2 is monotonic increasing in the (p, u) plane (p >^0). In addition, it is easily derived from (2.9) that 



lim u = —oo for the 2-shock curve S 2 , which implies that S 2 has the u-axis as its asymptotic line. It 

p—>0+ 

can also be derived from (2.8) that lim u = —00 for the 1-shock curve Si. 

p—>- + oo 

In the phase plane (p > 0, u € R), through point (p_, u_), we draw the elementary wave curves Ri, 
R2, Si and S2, respectively. Then the phase plane is divided into four regions I, II, III and IV(p_,u_) 
(see Fig. 1). 

By the analysis method in phase plane, for any given state (p+, u+), one can construct the Riemann 
solutions as follows: 

(1) {p+,u+) e I{p-,u-) : Ri +R2; 

( 2 ) (p-i-,u-i-) e //(p_,u_) : Ri + S2; 

( 3 ) (p+,M+) G Ill{p-,u-) : Si + R 2 ; 

(4) {p+,u+)GlV{p-,U-):Si + S2. 



u 


Fig. 1. Curves of elementary waves. 


Thus we have proved the following result 

Theorem 1. For Riemann problem (1.1), (1.2) and (2.1), there exists a unique entropy solution, which 
consists of shock waves, rarefaction waves, and constant states. 


3. Riemann problem for the transport equations 

The Riemann solutions to the transport equations in zero-pressure flow were presented by Sheng 
and Zhang in [14]. The Riemann problem to the transport equations are 

I Pt + (pu). = 0, 

I (pu)t -I- (pU^)a: = 0 


with initial data 


(p, u){x, 0) = (p±gU±), ±x > 0. 


(3.2) 



The system has a double eigenvalue 


X = u 


and only one right eigenvector 

1 ^ = (r, 0)^. 


By a direct calculation, 


VA-y = 0. 


Thus (3.1) is nonstrictly hyperbolic and A is linearly degenerate. 

As usual, we seek the self-similar solution 

[p,u){x,t) = {p,u){i), C=|- 

Then Riemann problem (3.1) and (3.2) is reduced to the following boundary value problem of ordinary 
differential equations: 

( 

-Cp? + (p-w)c = 0, 

-C(pm)? + (p«^)? = 0, 

with (p,m)(±oo) = (p±,m±). 

For any smooth solution, system (3.3) can be written as 



= 0 . 


(3.3) 


(3.4) 


(3.5) 


It provides either the general solution (constant state) 

(p, m)(C) = constant (p / 0) 

or the singular solution 

( 

p = 0, 

u = i, 

which is called the vacuum state (see [14]), where u(^) is an arbitrary smooth function. 

For a bounded discontinuity at ^ = ct, the Rankine-Hugoniot condition holds: 

-cr[pl -b [pu] = 0, 

-a[pu\ -b [pu^l = 0, 

where [g] — q+ — q- denotes the jump of q across the discontinuity. By solving (3.6), we obtain 

J : ^ = a = U-{= X-) = u+{= \+), (3.7) 


(3.6) 


which is a contact discontinuity. It is a slip line and just the characteristic of solutions on both its sides 
in (a:,t)-plane. 

The Riemann problem (3.1) and (3.2) can be solved by contact discontinuities, vacuum or delta 
shock wave connecting two constant states (u±,ii±I. 



For the case «_ < u+, there is no characteristic passing through the region U-t < x < u+t and the 
vacuum appears in this region. The solution can be expressed as 

{ {p-,U-), —00<X<U-, 

(0,^), w_ < ^ < u+, (3-8) 

(p+,M+), «+ < C < +00- 

For the case u- = u+, it is easy to see that the constant states {p±,u±) can be connected by a 
contact discontinuity. 



Fig. 2. Characteristics overlap domain. 


For the case U- > u+, the characteristic lines originating from the origin will overlap in a domain 
fl, as shown in Fig. 2. So, singularity must happen in fl. It is easy to know that the singularity is 
impossible to be a jump with finite amplitude because the Rankine-Hugoniot condition is not satisfied 
on the bounded jump. In other words, there is no solution which is piecewise smooth and bounded. 
Motivated by [14], we seek solutions with delta distribution at the jump. 

To do so, a two-dimensional weighted delta function w{s)Sl supported on a smooth curve L = 
{(t(s), a;(s)) : a < s < 6} is defined by 

{w{s)5l,(p) = f w{s)ip{t{s),x{s))ds (3.9) 

J a 

for any ip € C^{R x R+). 

Let us consider a solution of (3.1) and (3.2) of the form 

{ (p-,M_), x<crt, 

{w{t)S{x — at),a), X = at, (3.10) 

{p+,u+), x>at, 

where cr is a constant, w{t) G C^[0,-boo), and 5(-) is the standard Dirac measure. x{t), 'w{t) and a 
are the location, weight and velocity of the delta shock, respectively. Then the following generalized 



Rankine-Hugoniot conditions hold: 


_ _ 
dt ^ 

^ = a[p] - [pu], (3.11) 

= a[pu] - [pu% 

where [p] = p+ — p_, with initial data 

(x,«;)(0) = (0,0). (3.12) 

In addition, to guarantee uniqueness, the delta shock wave should satisfy the entropy condition: 

u+ < a < U-. 


Solving the system of simple ordinary differential equations (3.11) with initial data (3.12), we have, 
when p_ = P+, 

x{t) = 2 ^ {p-u- — p+u+)t, 

^ = ^(w- +«+); 


when p- ^ P+, 


(+\ y/P—^— “h \/P+^4- , I -/ \, 

x{t) = -^——— t, w(t) = VP-P+(m- - u+)t, 

_ y/p^U- + y/^U+ 


4. Formation of ^-shocks 

In this section, we study the formation of 5-shock waves in the Riemann solutions of system (1.1) 
and (1.2) in the case {p+,u+) € IV{p-,U-) with U- > u+ as both the pressure and the magnetic field 
vanish. 


4.1. Limit behavior of Riemann solutions as ki,k 2 —)■ 0 


When {p+,u+) £ IV{p-,U-), for each pair of fixed fci > 0 and k 2 > 0, suppose that (p*,-u*) is the 
intermediate state connected with (p-,u-) by a 1-shock with speed cri and (p+,u+) by a 2-shock 
S 2 with speed CT 2 . Then it follows 


Si : 


CTl = M_ — p, 


M» = M_ — 


/ 1 / ki 

' y p»p- p»p- 


k'^{pt+P-) 

2 II 




rr 

f fci I fc^(p.+p_)\ 

V 

yp,p- 2p y 


S2 : 




(4.1) 


(4.2) 


M+ = U* -I- . 





The addition of (4.1) and (4.2) gives 


U- — u+ = 


1 1 
P- P* 


,, 1 1 ,, mpi-pi) 

fcl(-) H-- 

P- p* 2^ 


1 


P+ P* 


+\/t— — ~) + 


1 


P+ P* 


2p 


p* > p±. 


(4.3) 


For any given p± > 0, if lim pt, = M £ [max(p_,p+),+oo), then by taking the limit ki,k 2 —>■ 0 in 

ki ,k2^0 

(4.3), we have U- — u+ — 0, which contradicts with m_ > m+. Therefore, lim p* = +cx3. Letting 

ki ,k2^0 

fci, fc 2 —>■ 0 in (4.3), we obtain the following result. 


Lemma 1. 


Lemma 2. 


klM-^0 {y/^ + y/p^)^ 


lim M« = hm CTi = hm (T 2 = - -=- — = a, 

fcl,/c2—^0 ki,k2—^0 kiyk2 —^0 






lim 

k\ ,k2—^0 


p,dx = (o-[p] — [pu])t = y/p-p+{u- — u+)t = w{t). 


(4.4) 

(4.5) 

(4.6) 


Proof. Letting fci, ^2 —>■ 0 in (4.1) and noting Lemma 4.1, we have 


1 - 1- /I 1 /i / 1 1 1 , ^2(P* -P-) 

lim u*=u-— lim A - \ ki( -)H--- 

fci,/c2^o ki,k2^o \ P- p* \ p- p* 2p 


_ ^ 1 p-p+{u- — U+)^ _ a/^U- + a/P+U+ _ ^ 


p-V + + 

From the first equation of (4.1), by Lemma 4.1, we obtain 


(4.7) 


lim ai = u- — lim i ^ + 
fcj,fc2^o fci,fc2^o y p_ 




2 p 


p+{u-—u+)^ y/p^U- + y/p+U+ 

= «_ — ,/ —-i=- — -==- - -- = cr. 


iV^ + + '/P+ 

From (4.2) and (4.7), we can easily get 


(4.8) 


1 Pi ^ 2 P+(~ + ~) 

lim CT 2 = lim m* + lim 1 / — H- = n. 


k\,k2—^0 ki,k2—^0 ki,k2—^0 V p* 


2 p 


(4.9) 


Thus it can be seen from (4.8) and (4.9) that when fci, P 2 —>■ 0, the two shocks Si and S 2 will coincide 
whose velocities are identical with that of the delta shock wave of the transport equations with the 
same Riemann initial data (p±,M±). 

Using the Rankine-Hugoniot conditions (2.7) for Si and S 2 , we have 


CTl(p* — p_) = p*M* — P-U-, 
(T2(p+ — p*) = P+U+ — p,u^. 


(4.10) 


Then from (4.8) and (4.9) it follows that 

lim ((Ti — o- 2 )p* = lim (p+M+— p_u_ + (Jip_ — cr 2 p+) = [pu] — (t[p]. 

fcl,fc 2->0 fcl,fe 2-»0 Q 


(4.11) 



This implies that 


lim / p,dx — {a[p] — [pu])t = yfpZp^{u-— u+)t = w{t). (4-12) 

The proof is completed. 

Remark 1. From the above results, it can be seen that the limit of the Riemann solution of system 

(1.1) and (1.2) as k\,k 2 —>■ 0 in the case (p+,M+) G IV{p-,U-) is just the delta shock solution of 

(3.1) -(3.2) when U- > u+. 

4.2. (5-shocks and concentration 


Now, we give the following results which give a very nice depiction of the limit in the case (p+,«+) G 
IV{p-,u-). 

Theorem 2. Let U- > m+ and {p+,u+) G IV{p- ,U-). For any fixed ki,k2 > 0, assuming that {p,u) 
is a solution containing two shocks Si and S 2 of (1.1)-(1.2) with Riemann initial data (2.1), constructed 
in Section 2, it is obtained that as ki,k 2 —>■ 0, {p,u) converges in the sense of distributions, and the 
limit functions p and pu are the sums of a step function and a (5-measure with weights 

{<j[p] — [pu])t and {(j{pu] — [pu^])t, 


respectively, which form a delta shock wave solution of (3.1) with the same Riemann initial data 

(p±, Ui). 

Proof. Let ^ = x/t. Then for any fixed fci > 0 and k 2 > 0, the Riemann solution to the isentropic 
magnetogasdynamics equations for Chaplygin gas (1.1)-(1.2) can be written as 


(p,m)(C) = 


(p»,M*), (Tl < (f < (T2, 
{P+,U+), C > 0-2, 


(4.13) 


which satisfies the following weak formulations: 






Piomd^ = 0 


(4.14) 


and 

J {^-ui0)p{0ui0'>P'i0d^- J )V’'(C)ciC+ J P(C)M(?)^(6)rfC = 0 (4.15) 

for any test function tp G C“(—00, -|-oo). 

The first integral on the left-hand side of (4.15) can be decomposed into 

( PCTl p(T2 p + OO 


-f 


-f 




(4.16) 



The sum of the first and the last terms in (4.16) is 


^ + oo 




= p_M_CTi'i/)(cri) — P-U^_'lp{a\) — P-U- 


-p+U+a2'lp{o2) + p+u\'lp{(T2) - P+U+ 


> —oo 
0 + 00 




Taking the limit fci, —>■ 0 in (4.17) leads to 


lim 

ki ,^2—^0 


r T r 

' —oo j (T 2 


(C - «(C))p(0«(C)V''(0rfC 


(4.17) 


= ([pw^] - +pu]++) - / (PoWo)(C - + ■ 

^ —OO 

where (poUo)(C) = P-U- + [pu]H{^) and H is the Heaviside function. 
For the second term in (4.16), integrating by parts again, we obtain 

/ (C -'«(C))p(C)«(C)V''(6)rfC = / (C - w*)p*M*^/''(0rfC 


(4.18) 


= -pM,{pp{o 2 ) - i’+l)) + P*u,{a 2 tp{a 2 ) - criV'(cri)) - p*M* 


= —M*P»((T2 — (Tl) 


V>(cr2 ) - ■!/>(cti) (T 2 'lp{o- 2 ) - cri'lp{ai) 

-W*-h 


1 /■' 


(72 — (71 


(72 — (71 


(72 — (71 




(4.19) 


Taking the limit ki,k 2 —>■ 0 in (4.19), noting (4.11) and the fact that both tp € C“(—oo,+oo) and 
lim M* = lim ai = lim (72 = cr, we deduce that 

/ci,/c2—>-0 ki,k2^0 ki,k2—^0 


lim [ {^-u{^))p{^)u{^)+i^)d(,^a{[pu]-a[p]){aip'{a)-aip'{a)-ip{a) + p!{a))^0. (4.20) 

Similarly, the first integral on the left-hand side of (4.15) can be decomposed into three parts as 

ki , k^piOf 




(4.21) 


which equals to 




"" “ V'(o-i)) - ■^^(V’(f^2) - ^/'((7i)) - - -|Y^)^('^2). (4.22) 


fcip' 


fci fcip+' 


Taking the limit ki,k 2 —>■ 0 in (4.22), by Lemmas 4.1-4.2, we have 


lim 

All ,A:2~^0 

Summarizing (4.18), (4.20) and (4.23) leads to 

/* + oo 



kliPiOf) 

/-oo Vp(C) 

2p J 


= 0 . 


lim 

/ci ,fc2—>-0 


((pM)(C) - (poWo)(C - + + iOd^ = +[pu] - [pu ])ij{(r) 
11 


(4.23) 


(4.24) 



which is true for any ip € C'“(—oo, +oo). 

As done previously, we can obtain the limit for the first integral on the left-hand side of (4.14) as 

/ + 00 p(T ^ + 00 

(C -'«(C))p(C)^'(C)c*? = (M - cr[p])^(cf) - / 

/ / 

•OO j —OO j <7 


= i[pu] - cr[p])V'(o-) - 


po(C- 


where po(C) = P- + [p]^(C)- Then returning to the formulation (4.14), we have 

P + OO 


lim 

ki ,k2 —^0 


(p(C) - Po(C - cr))^(6rfC = {^[p] - [Pu])'ip{(j), 


(4.25) 


(4.26) 


which is true for any ip G C”)—oo, -too). 

Finally, we study the limits of p and pu as fci, ^2 —>■ 0, by tracing the time-dependence of weights 
of the (5-measure. Let <p{x,t) G oo,-l-oo) x [0,-|-oo)), then we have 


^ + oo ^ + oo 


r + oo 


lim 

fci ,^2 —>-0 


p{x/t)(p{x,t)dxdt = lim 

ki,k2—^0 


^+oo 


]dt. 


Regarding t as a parameter and applying (4.26), one can easily see that 


lim 

k\ ,k2—^0 


p-\-oo p + oo 

/ P{^)H^t,i)d^= Po{^-(^)cp{^t,t)d^ +{a[p]-[pu])(p{at,t) 

' —OO J —OO 


P+OO 


Po (I - <P{x, t)dx + {a[p] - [pu])(p{Gt, t), 

Substituting (4.28) into (4.27) and noting po(f ~ ~ Po(® — <xt), we have 


Hoo /*+oo 


1-00 /*-|-oo 


lim 

ki ,k2—^0 


p{x/t)(p{x,t)dxdt = 


po{x — at)(p{x,t)dxdt 


(4.27) 


(4.28) 


+ / t{a[p] — [pu])(p{at,t)dt. (4.29) 

Jo 

By definition (3.9), the last term on the right-hand side of (4.29) equals to {wi{t)5s, <p{-,-)): where 

wi(t) = {g[p] - [pu])t. 

With the same reason as before, we arrive at 


hoo /»-|-oo 


|-oo /»-|-oo 


lim 

1^2 


p{x/t)u{x/t)(p{x,t)dxdt = 


{poUo){x — at)(p{x,t)dxdt 


+ / t{a[pu] — [pu‘^])(p{o't,t)dt. 

Jo 

The last term on the right-hand side of (4.30) equals to {w 2 {t)Ss, (p{-, •))) where 

W2{t) = {a[pu] - [pu^])t. 


(4.30) 


The proof is completed. 


5. Formation of vacuum st^^es 



In this section, we study the formation of vacuum states in the Riemann solutions of system (1.1) 
and (1.2) in the case (p+,M+) € 7(p_,u_) with u_ < u+ and p± > 0 as both the pressure and 
the magnetic field vanish. In this case, we know that the Riemann solution consists of a backward 
rarefaction wave R\, a forward rarefaction wave R 2 and an intermediate state besides two 

constant states (p±,M±), which are as follows 





For any given p± > 0, if lim p, = K £ (0, min(p_, p+)], then by 
fcl , fc 2—^0 


(5.3) 



it follows that 



Therefore, by the squeeze theorem in multivariable calculus, we arrive at 



Similarly, we can obtain that 



Combining (5.3), (5.6) and (5.7), we have U- — u+ = 0, which contradicts with U- < u+. Therefore, 
lim p* = 0, which implies that a vacuum occurs. From (5.1), one can see that 

ki ,k2—^0 
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< U- — 


ki k^P 

4 + — + 


It can be derived from (5.4) that 




P*<P<P- 


(5.8) 


U- - \ — + 


fci 

p^ 


kip 


+ 


fci I 

TT H-T" 


-ds 




p2 p 


^ _i_ 


+ \/ -J + 




fcip- 


In I \ / — + 


fci 


+ 


fc|p_ 


klp- 


Inp 


^1 , fciP- , klp- 


P- 


P 


P 


■In 


I k\ klp- 

+ 

P- P 


fcjp- 

P 


fc|p_ 


lnp_. 


(5.9) 


The uniform boundedness of p(^) with respect to ki, k2 in this case leads to 


lim 

k\ ,^2—^0 


fei fcjp 

9 ~r I 

p^ p 


r- v'4+ 


-ds I = lim 

k\,k2 —^0 


9 

p2 p 


= u-. (5.10) 


Then, by the squeeze theorem in multivariable calculus, we have lim Ai — U-. Similarly, we can 

/ci ,fc2—>-0 

obtain that 


lim A2 = u+ and lim u(^) = for ^ £ (m_,m+). 

fci , fc 2—^0 ^ 15 ^ 2—^0 


(5.11) 


Then from above we have proved the following results. 


Theorem 3. In the case (p+,M+) £ /(p_,M_) with u_ < u+, as fci,fc 2 — >■ 0, the vacuum state 
occurs and two rarefaction waves Ri and R 2 become two contact discontinuities u = U- and u = u+, 
respectively, connecting the constant states (p±,u±) with the vacuum (p = 0). 

Theorem 4. In the case (p+,u+) £ 7(p_,M_) with U- < u+, as ki,k 2 — >■ 0, the limit of the Riemann 
solution of (1.1) and (1.2) with initial data (2.1) is just the Riemann solution of the transport equations 
(3.1) for zero pressure flow with the same initial data, which contains two contact discontinuities 
xjt = u± and a vacuum state besides two constant states. 
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